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ABSTRACT 

By using a one-zone chemical and spectrophotometric evolution model of a 
disk galaxy undergoing a dusty starburst, we investigate, numerically, the optical 
spectroscopic properties in order to explore galaxy evolution in distant clusters. 
We adopt an assumption that the degree of dust extinction (represented by Ay) 
depends on the ages of starburst populations in such a way that younger stars 
have larger Ay (originally referred to as selective dust extinction by Poggianti & 
Wu 2000). In particular, we investigate how the time evolution of the equivalent 
widths of [OII]A3727 and H5 is controlled by the adopted age dependence. This 
leads to three main results: (1) If a young stellar population (with the age of 
~ 10^ yr) is more heavily obscured by dust than an old one (> 10* yr), the 
galaxy can show an "e(a)" spectrum characterized by strong US absorption 
and relatively modest [Oil] emission. (2) A dusty starburst galaxy with an e(a) 
spectrum can evolve into a poststarburst galaxy with an a-|-k (or k-|-a) spectrum 

0. 2 Gyr after the starburst and then into a passive one with a k-type spectrum 
1 Gyr after the starburst. This result clearly demonstrates an evolutionary link 
between galaxies with different spectral classes (i.e., e(b), e(a), a-|-k, k-|-a, and 
k). (3) A dusty starburst galaxy can show an a-|-k or k-|-a spectrum even in the 
dusty starburst phase if the age- dependence of dust extinction is rather weak, 

1. e., if young starburst populations with different ages (< 10^ yr) are uniformly 
obscured by dust. 

Subject headings: galaxies: clusters - galaxies: formation - galaxies: ISM - 
galaxies: infrared - galaxies: interaction - galaxies: structure 
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1. Introduction 

Spectroscopic properties of galaxies have provided valuable information on different 
aspects of galaxy evolution, such as star formation histories along the Hubble sequence, 
chemical evolution of gas and stars, and the physical conditions associated with AGN 
and starburst activity in the central regions of galaxies (e.g., Kennicutt 1998, 1999). In 
particular, spectroscopic studies of galaxies at optical wavelengths have played vital roles 
not only in determining the current and the past star formation rate but also in deducing 
the physical mechanisms relevant to the time evolution of star formation (see Kennicutt 
1998 for a recent review). The change in the galactic global star formation rate with redshift 
has recently been investigated for a variety of different distant star-forming objects by 
many authors (e.g., Madau 1996; Steidel et al. 1999; Ellis 1997). These recent works have 
suggested that optical data alone cannot provide a comprehensive view of the star formation 
properties of distant galaxies owing to heavy dust extinction. A major and important 
challenge is to properly assess the extent to which the stellar hght is affected by dust 
extinction. Future near-infrared spectroscopic studies of the redshifted Ha emission from 
distant, dust-obscured galaxies is crucial to revealing the degree of extinction and reddening 
and thereby clarifying the real star formation histories of these galaxies (Kennicutt 1999). 

Detailed spectroscopic studies of distant rich clusters (Dressier & Gunn 1983, 1992; 
Couch & Sharpies 1987; Poggianti et al. 1999) have provided important clues as to the 
nature and evolution of galaxies in these dense environments. In particular. Dressier & 
Gunn (1983, 1992) discovered galaxies with no detectable emission lines but with very 
strong Balmer line absorption (the so-called "E+A" galaxies) and concluded that such 
galaxies were "poststarburst" objects seen within ~ 1 Gyr of a severely curtailed burst of 
star formation. Furthermore, Couch & Sharpies (1987) made the first attempts to establish 
an evolutionary link between star-forming, poststarburst, and passively evolving galaxies 
observed in distant clusters. 

These attempts have been progressed significantly further by Dressier et al. (1999) 
and Poggianti et al. (1999) through their detailed modeling of an extensive data base of 
spectroscopic observations they obtained in 10 rich clusters in the range 0.37 < z < 0.56. 
The wealth of data allowed them to delineate the different spectral classes pertinent to 
distant cluster populations much more clearly, with three distinct classes of emission line 
object [e(a), e(b), e(c)] being identified on the basis of the strength of [OII]A3727 emission 
and US absorption, two poststarburst classes (a+k and k+a) identified amongst galaxies 
with no detectable [OH] emission but enhanced H5 absorption, and a passive "k" class 
where neither of these features had any prominence. Attempts were made to identify the 
evolutionary pathways between these different classes of object in order to identify the 
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most likely progenitor to the sizeable population of poststarburst galaxies. Galaxies with 
e(a) spectra characterized by strong Balmer absorption [EW(H5) > 4A] and relatively 
modest [Oil] emission, were previously interpreted (in the context of dust-free models) as 
poststarburst galaxies with some residual star formation (Couch et al. 1998). However, 
Poggianti et al's detailed spectral analysis and comparison with Liu & Kennicutt's (1995) 
sample of spectra for present-day merging galaxies suggested that the e(a) types were more 
hkely dusty starburst galaxies. This suggestion has subsequently been reinforced by the 
observation that a significant number of e(a) galaxies are either luminous infrared galaxies 
(LIRGs) or ultra-luminous infrared galaxies (ULIRGs), objects which are often dominated 
by dusty starbursts (Poggianti & Wu 2000). 

Independent evidence for the presence and importance of dust in determining the 
observed properties of distant cluster galaxies has been presented by other authors. Small 
et al. (1999) found that 5 out of 10 galaxies classified as either a-|-k or k-|-a types by 
Dressier et al. (1999) in Abell 851 (at z = 0.42) were detectable in deep 1.4 GHz continuum 
maps obtained with the VLA. With this radio emission most likely attributable to starburst 
activity, the lack of any such signature in the optical spectrum could only imply that it 
must be totally obscured by dust. Further to this, Owen et al. (1999) investigated the radio 
properties of galaxies in a rich cluster at z ~ 0.25 (A2125) which had large blue galaxy 
fractions (0.19) and found that the optical line luminosities (e.g., Hq;-|-[NII]) were often 
weaker than one would expect for the star formation rates implied by the radio emission. 
They accordingly suggested that dust obscuration, which is larger than is usually found 
locally, hides most of optical hght from the star-forming regions. 

Although these observational results have identified the possible importance of dust 
and dusty starbursts in distant cluster populations, little theoretical work has been done to 
address the formation and evolution of such galaxies. The purpose of this paper is to report 
the results of a numerical investigation of the spectroscopic evolution of dusty starburst 
galaxies and possible evolutionary links between galaxies of different spectral class. Using a 
one-zone chemical and spectrophotometric evolution model, we consider in particular, the 
degree to which dust extinction (represented by Ay) depends on the age of the starburst 
populations and thereby try to clarify the time evolution in the equivalent width of the 
H(5 absorption - and [OH] emission - line features. The main questions that wc address 
are: (1) When and how does a disk galaxy hosting a dusty starburst develop an e(a)-type 
spectrum? (2) What are the requisite physical conditions for the formation of e(a) galaxies? 
(3) Can a dusty galaxy show an a-|-k or k-|-a spectrum even in its starburst phase? (4) Are 
there any evolutionary hnks between the e(a), e(b), e(c), a-|-k, k-|-a, and k spectral types? 

A key feature of our approach is to adopt the selective dust extinction framework 
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originally proposed by Poggianti & Wu (2000) whereby the extinction suffered by a stellar 
population is inversely proportional to its age with the youngest popiilations having 
the greatest extinction. We confirm that such age-dependent extinction is essential to 
reproducing the spectral properties of the putative dusty starburst population in distant 
clusters. 

The layout of this paper is as follows: In §2, we summarize the numerical models 
used in the present study and describe briefiy the methods for deriving spectral energy 
distributions (SEDs) affected by internal dust extinction. In §3, we present our numerical 
results on the time evolution of dusty starburst galaxies in the EW([011])-EW(H(5) plane. 
In §4, we discuss the possible evolutionary connections between the different spectral types 
(e.g., e(a), e(b), and k+a) in the context of our dusty starburst models. The conclusions of 
the present study are given in §5. In the present paper, each spectral class is defined as 
follows: e(a) for galaxies with 5 A< EW([OII]) < 40 Aand 4 A< EW(H(5), e(b) for 40 A< 
EW([OII]) and EW(H(5) < 4 A, e(c) for 5 A< EW([OII]) < 40 Aand EW(H5) < 4 A, a+k 
for EW([OII]) < 5 Aand 7 A< EW(H5), k+a for EW([OII]) < 5 Aand 3 A< EW(H(5) < 7 
A, k for EW([OII]) < 5 Aand EW(H5) < 3 A. 

2. Model 

We adopt a one-zone chemical and spectrophotometric evolution model of a disk galaxy 
with a starburst, and thereby investigate the time evolution of the line strength of 115 
absorption and that of [Oil] emission, both of which are important properties for classifying 
the spectral types of galaxies in distant clusters (Dressier et al. 1999). The present model is 
essentially the same as that adopted in our previous studies (Shioya & Bekki 1998, 2000), 
with the details of our treatment of chemical evolution given in Shioya & Bekki (1998). 
Accordingly, we describe only briefiy the present spectrophotometric model. 

One of the main differences between our model and previous ones (e.g., Couch & 
Sharpies 1987; Newberry et al. 1990; Chariot & Silk 1994; Jablonka & AUoin 1995; 
Poggianti & Barbaro 1996; Poggianti et al. 1999) is that we explicitly include the effects of 
dust on both the photometric and spectroscopic properties of stellar populations as also the 
gaseous emissions. With this new development to our code, we can clearly demonstrate how 
the dust extinction affects the time evolution of the fine strengths of R6 and [Oil] (EW(H5) 
and EW([OII])) in dusty starburst galaxies. 

We follow the chemical evolution of galaxies by using the model described in Matteucci 
& Tornambe (1987) which includes metal-enrichment processes of Type la and II supernovae 
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(SNIa and SNII). We adopt the Salpeter initial mass function (IMF), (/)(m) oc m~^-^^, with 
upper mass hmit Mup = 120 M^. and lower mass limit M\o„ = O.IMq. We calculate the 
photometric properties of galaxies as follows: The monochromatic flux of a galaxy with age 
T, Fx{T), is described as: 

Fa(T) = r FsspAZ,T-t)i/j{t)dt, (1) 
Jo 

where Fssp,\{Z,T — t) is the monochromatic fiux of a single stellar population with age 
T t and metallicity Z, and ipit) is the time-dependent star formation rate described 
below. In the present study, we use the spectral library GISSEL96 which is the most recent 
version of that published by Bruzual & Chariot (1993). 

The star formation history of a disk galaxy is characterized by three epochs. The first 
is the epoch of galaxy formation (Tform) during which the interstellar gas is consumed by 
star formation at a moderate rate. In the present study, Tform is for all models. The 
second is Tgb at which the starburst begins in the disk. Here is considered to be a free 
parameter. The third is Tend at which all star formation ceases is defined as the epoch at 
which the stellar mass fraction reaches a value /final in our models; /gnai is also assumed to 
be a free parameter which controls the duration of starburst. Throughout the evolution of 
disk galaxies, the star formation rate is assumed to be proportional to the gas mass fraction 

m = ku (2) 

where k is a parameter which controls the star formation rate. This parameter k is given as 
follows: 

/cdisk for < T < Tsb, 
k^l /csb for T,b < T < Tend, (3) 

, for Tend < T. 

In the following, we refer to these three different epochs as the "prestarburst" , "starburst" , 
and "poststarburst" phases, respectively. 

The value of the /c^isk parameter, which controls the star formation histories of disk 
galaxies, is suggested to vary for different Hubble types (e.g., Arimoto, Yoshii, & Takahara 
1992). For example, fcdisk = 0.225 Gyr~^ corresponds to a plausible star formation rate for 
Sb disks (e.g., Arimoto et al. 1992). In comparison, observations of starburst galaxies (e.g., 
Planesas et al. 1997) would suggest that ksb is an order of magnitude larger than /cdisk (i-e., 
ksh — 2.25 Gyr~^). By changing the value of kg^, we can investigate how the strength of a 
starburst influences the evolution in spectral type within a dusty galaxy. 

We adopt the selective dust extinction model investigated in detail by Shioya & Bekki 
(2000). In this model, the amount of dust extinction associated with the starburst is 
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assumed to depend on the age of stellar population: 



Av = Av,o exp{(T - t)/T} for T^b < T < Tend, 



(4) 



where Ay^ and r are free parameters which control the the extinction at any given time 
during starburst. Based on the monochromatic flux derived in equation (1) and the value 
of Ay given by equation (4), we calculated a 'dust-absorbed' SED for galaxies of different 
age T using the extinction law inferred by Calzctti ct al. (1995) and adopting the so-called 
screen model. To derive fluxes for various gaseous emission lines (B.6 and [Oil]), we flrst 
calculate the number of Lyman continuum photons, A^Ly, by using a non-absorbed SED. 
If all the Lyman continuum photons are consumed in ionizing the surrounding gas, the 
luminosity of H/3 can be calculated according to the relation: 



(Leitherer & Heckman 1995). To calculate luminosities of other emission lines, e.g., [Oil] 
and H5, we use the luminosity ratios relative to H/9 tabulated in PEGASE (Fioc & 
Rocca-Volmcrangc 1997), which arc calculated for an electron temperature of 10,000 K and 
an electron density of 1 cm~^. The strength of absorption lines arc derived using GISSEL96 
(Bruzual & Chariot 1993). Thus the SEDs derived here comprise the following components: 
stellar continuum, gaseous emission, and stellar absorption. 



Shioya & Bekki (2000) have already demonstrated the importance of selective dust 
extinction in reproducing c(a) type galaxy spectra by comparing the model results with and 
without dust extinction. Furthermore, they pointed out that a dusty starburst galaxy with 
selective dust extinction evolves from an e(b), to an e(a), to an a-|-k, to a k-|-a, and finally 
to a k type galaxy. Since they have already indicated the fundamental role of the selective 
dust extinction in spectroscopic evolution of dusty starburst galaxies, we detail here the 
parameter dependence of the spectroscopic evolution in the EW([0II])-EW(H(5) plane. To 
demonstrate the importance of each of the six free parameters, Ayfl, r, Tgb, k^isk, ksh, and 
/final, we firstly compute a fiducial model (Model 1) with Avfi = 5.0 mag, r = 1.0 x 10® yr, 
= 7.64 Gyr, A;disk = 0.225 Gjr-\ k.^ = 2.25 Gyr'^, and /final = 0.95. Model 1 can then 
be compared with a model without dust extinction (Model 2), to show the typical behavior 
of selective dust extinction in the evolution of starburst galaxies on the EW([0II])-EW(H5) 
plane. The dependence of the evolution on the six parameters, Avfi, t, ^sb, ^disk, ^sb, and 
/final is then investigated via the computation of 17 models with a range of parameter 



L(H/5)(erg s"^) = 4.76 x 10~^^iVLy(s"^) 



(5) 



2.1. Main points of analysis 



values: 
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Model 1, 3, 4, and 5 with different Ay^, 

Model 1,6, and 7 with different r. 

Model 1, 8, 9, 10, and 11 with different Tgb, 

Model 1, 12, and 13 with different fcdisk, 

Model 1, 14, and 15 with different ksh, 

and Model 1, 16, and 17 with different /cfinai- 

The parameter values are given in Table 1. We assume that the initial gas mass of the 
model galaxy is IO^^Mq. 

3. Results 

3.1. The fiducial model 

Figure 1 shows the star formation history and the color evolution (in B — V) for 
the fiducial model (Model 1) with fcdisk = 0.225 Gyr~^ (corresponding to Sb-type star 
formation history). During the starburst, the B — V color becomes rather blue (0.2 mag 
at T = 8.56 Gyr) due to the formation of massive stars. However, the color is not so 
blue compared with the model without selective dust extinction (Model 2) because of the 
heavy dust obscuration suffered by the young stellar population in the fiducial model. This 
result implies that selective dust extinction affects the global colors as well as the spectral 
characteristics of starburst galaxies. 

Figure 2 shows the evolution of the emission and absorption lines by plotting their 
strength at T = 1, 7.64, 8.56, 8.58, and 13.0 Gyr for the fiducial model. Before the starburst 
{T < 7.64 Gyr), the model clearly shows a number of strong emission lines, whereas the 
absorption lines (e.g., the B.S absorption line) are rather weak as a result of the continuous 
and moderate star formation. During the starburst (7.64 < T < 8.58 Gyr), the emission 
lines do not become so remarkably strong; nevertheless, the star formation rate becomes 
more than a factor of five larger than that of the prestarburst phase. This is because the 
ionizing photons from very young OB stars are very heavily obscured by the dust. After the 
starburst (T > 8.58 Gyr), the emission fines become even less remarkable due to the death 
of the young, massive stars and the intermediate-aged ones as well. On the other hand, 
the stellar absorption fines (e.g., R6) become prominent after the starburst {T > 8.58 Gyr) 
owing to the presence of an older stellar population. 
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Figure 3 describes the time evolution of EW([OII]) and EW(H5), which clearly 
demonstrates the importance of the selective dust extinction in the spectroscopic evolution 
of dusty starburst galaxies. The equivalent width of [Oil] does not increase and is well 
below 40 A during the starburst. The equivalent width of H5 becomes rapidly smaller 
just after the starburst begins, although the H5 line can still be observed as an absorption 
line. A few xlO^ years after the starburst begins, the equivalent width of B.6 becomes 
larger than 4 A. Consequently, this model reproduces the properties of the e(a) class with 
EW(H5) > 4 Aand relatively modest [Oil] emission, EW([OII])> 5 A; the location in the 
EW([0II])-EW(H5) plane is shown in Figure 4. The physical reason for the successful 
reproduction of the model is clearly described in our previous paper (Shioya & Bekki 2000). 
The time scale during which the model is located within the e(a) region is estimated to be 
~ 0.7 Gyr. These results thus confirm the earlier suggestion of Poggianti & Wu (2000) of 
the vital role of selective dust extinction in reproducing e(a) spectra. 

3.2. Pcirameter dependences 

Figure 5 shows the time evolution of spectroscopic properties in the EW([0II])-EW(H5) 
plane for models with different values of Ayfl (1 < Ay,o < 10-0 mag for Model 2, 3, 4, 
and 5). Since parameters other than Ay^ are fixed at the same values adopted for the 
fiducial model, this figure demonstrates how the evolution of dusty starburst galaxies in the 
EW([011])-EW(H5) plane differs according to the amount of dust extinction. As can be 
seen in the figure, the locus on the plane is appreciably different between the models, which 
implies that the degree of dust extinction (represented by Ayfl) is an important influence 
on the observed spectroscopic evolution. During the starburst, the model with the larger 
Avfi shows weaker [Oil] emission, weaker R5 emission, and stronger R6 absorption. Thus 
the model with the largest Ayfi (e.g.. Model 5 with Ayfi — 10 mag) is located towards the 
bottom-right of the e(a) region during the dusty starburst. 

Figure 6 describes the evolution in the EW([0II])-EW(H5) plane for the models with 
different values of r (1 x 10^ < r < 1 x lO^yr for Model 1,6, and 7). Parameters other than 
r are fixed at the same values of the fiducial model. This figure clearly demonstrates that 
the time-scale during which young massive stars are heavily obscured by dust also governs 
the spectroscopic evolution. For the parameter range adopted here, a larger number of 
young massive stars with different ages are obscured by dust for the model with larger r, 
i.e., the age-dependent dust extinction is weaker and thus young stars with different ages 
(< 10^ yr) are more uniformly obscured by dust. Therefore, the model with larger r shows 
lower EW([OII]) and higher EW(H5) and, as a result, is located in the bottom- right hand 
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region of the EW([0II])-EW(H(5) plane during the starburst. We stress here that for the 
model with larger r (= 3 x 10^ yr and 1 x 10^ yr; Model 6 and 7, respectively), the time-scale 
during which the galaxies show an e(a)-type spectrum is appreciably shorter compared to 
the fiducial model. Furthermore, the galaxies show a k+a/a+k spectral signature in the late 
phase of the starburst. This result suggests that not only the galaxies with e(a) spectra but 
also those with k+a/a+k spectra arc possible dusty starburst candidates in distant clusters. 

For the models with r > 1 x 10^ yr, not only young massive stars (e.g., OB stars) - 
which emit most of the ionizing photons - but also late B and A stars are strongly affected 
by dust extinction. Hence a galaxy with a r value of this size docs not show strong B.6 
absorption during the starburst and poststarburst phases. Therefore the galaxy can be 
located in k+a region even in the dusty starburst phase. This can be contrasted with models 
with T < 1 X 10^ yr, where the emission lines from ionized gas regions around massive stars 
are not obscured so heavily by dust. As a result of this, even if we adopt a rather large 
value for Ayfi (i-e., more than 10 mag), neither the [Oil] nor B.5 lines are weakened. In this 
dusty starburst galaxy shows an e(b) spectrum during its starburst phase. 

These results on the r-dependences within our models confirms that age-dependent 
dust extinction is a crucially important factor for determining the spectroscopic properties 
of dusty starburst galaxies. The present onc-zonc model, however, cannot determine with 
any great certainty what controls the age- dependence (i.e., the parameter value of r) owing 
to the limitations of the model. We further discuss this problem in a later section. 

Figure 7 describes how the star formation history and the evolution of spectroscopic 
properties on the EW([0II])-EW(H5) plane are controlled by the parameter Tgb in dusty 
starburst galaxies. For models with a smaller Tgb, the starburst occurs at earlier epochs 
and the magnitude of the starburst is larger. This is because there is more gas to fuel 
starbursts at earlier times. The evolution of galaxies on the EW([011])-EW(H5) plane 
is not so remarkably different between the models with different Tgh (=4.46, 6.93, 7.64, 
8.47, and 9.45 Gyr: corresponding to Model 8, 9, 1, 10, and 11, respectively). The model 
with smaller Tgb (e.g.. Model 9), which represents a higher redshift starburst galaxy, shows 
discernibly lower EW([OII]) and higher EW(H(5) and thus is located at the bottom- right of 
the e(a) region during the starburst. This is primarily because a larger number of young 
massive stars are created and then obscured preferentially in the model with smaller Tg^. 
These results firstly imply that irrespective of the epoch at which a dusty starburst occurs, 
disk galaxies can exhibit an e(a) spectrum during this phase if the extinction parameters 
are within certain reasonable ranges. Secondly, they imply that the higher redshift dusty 
starburst galaxies are more hkely to be located to the right and to the bottom of the e(a) 
region in the EW([0II])-EW(H5) plane. It would be very interesting if future spectroscopic 
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observations could confirm this predicted redshift-dependent behavior in dusty starburst 
galaxies. 

Figure 8 describes how the evolution of starburst galaxies in the EW([0II])-EW(H5) 
plane depends on their prior star formation histories. The three models with Sa (Model 
13), Sb (Model 1), and Sc (Model 12) type star formation histories all pass through the 
e(a) region during the starburst, which implies that a dusty starburst galaxy can show an 
e(a) spectrum irrespective of its longer-term activity in its quiescent state. 

Figure 9 shows the dependence of the spectroscopic evolution on the parameter ksh 
controlling the strengths of starbursts in Sb type disk galaxies. Since parameters other than 
fcsb are the same between these models (with fcsb = 1-0, 2.25, and 5.0 Gyr~^, corresponding 
to Models 14, 1, and 15, respectively), this figure describes how the strength of a starburst 
in a galaxy determines the location of the dusty galaxy in the EW([0II])-EW(H5) plane. 
For the adopted parameter range (1 < ksh < 5.0 Gyr"^), the models show e(a) spectra 
during starbursts, although the location of each model in the e(a) region is different between 
these models: The models with larger kg^ (i.e., those with a stronger starburst) show 
stronger EW(H5) and are thus located to the right of the e(a) region. The two models with 
ksb — 2.25 and 5 Gyr~^ evolve from e(a) to a+k galaxies whereas the model with ksh — 1-0, 
corresponding to a galaxy with a rather weak and longer starburst, evolve from e(a) to 
k+a without passing through the a+k region. This result implies that dusty starburst 
galaxies that have undergone a rather weak starburst, are more likely to exhibit a k+a 
spectrum straight after the starburst. Clearly the strength of the starburst is important to 
the evolution of a dusty starburst in the EW([0II])-EW(H5) plane. 

Figure 10 demonstrates how the evolution is dependent on the duration of the 
starburst. There is about 1.8 Gyr difference in the starburst duration between the models, 
with /final = 0.90 (Model 16), 0.95 (Model 1, fiducial model), and 0.99 (Model 17). All 
three of the models pass through the e(a) region during the starburst, albeit at different 
locations. The models with smaller /gnai values (Model 1 and 14) - corresponding to those 
with the shorter starburst durations - evolve from e(a) into a+k and then into k+a types 
during the starburst and poststarburst phases. In contrast, the model with the largest /final 
value and thus the longest starburst duration, evolves directly from an e(a) to a k+a type 
without passing through the a+k spectral phase. (This dependence on /final is strikingly 
similar to that on ksh-) This result implies that a dusty starburst galaxy with the longer 
starburst duration is more likely to show the k+a spectra in the poststarburst phase. 

To summarize the above, weak [O II] equivalent width during the starburst is favoured 
by large dust obscuration, obscuration which is slow to clear relative to the characteristic age 
of the starburst. Furthermore, the location and movement within the EW([0II])-EW(H5) 
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plane of a dusty starburst galaxy during both the starburst and poststarburst phases has 
been demonstrated to be dependent on a number of parameters: Ayfl, r, T^b, kg},, A'disk- 
and /final- In turn, these parameters could depend on the relative spatial distribution of 
gas and stars, the radial variation in the metallicity and density of the gas, the dynamical 
evolution of gas and stars, and the effectiveness of supernovae feedback in the evolution of 
the interstellar medium for a dusty starburst galaxy. The one-zone model adopted here 
is, by its very nature, unable to determine the extent to which these important factors 
(e.g., the relative spatial distribution of gas and stars) are closely associated with the origin 
of the six parameters. This means that although the present study can indicate which of 
those parameters are important to the spectroscopic evolution of dusty starburst galaxies, 
it cannot so clearly reveal the origin of the spectroscopic properties of the galaxies. Thus 
future studies, in which the above five parameters are not treated as free parameters but as 
ones that can be predicted from a certain theoretical calculation arc required to understand 
more clearly the origin of each spectral type (e.g., e(a)) in the context of distant cluster 
galaxy evolution. 

4. Discussion 

4.1. Origin of selective dust extinction 

Calzetti et al. (1994) and Mayya & Prabhu (1996). have already proposed the idea 
that the older populations can be located at the outer part of molecular clouds and thus do 
not suffer so severely from dust extinction. Poggianti & Wu (2000) furthermore pointed 
out that stellar populations with the ages between a few 10^ and 1.5 xlO^yrs, which are 
responsible for the strong Balmer line absorption, have had sufficient time to drift away 
from the dusty molecular clouds where they were born. This selective obscuration of young 
stellar populations probably works effectively wherever a starburst occurs within molecular 
clouds. Both Shioya & Bckki (2000) and the present study have demonstrated that if 
younger starburst populations arc more heavily obscured by dusty gas than old ones in a 
galaxy, the galaxy shows an e(a)-type spectrum. This was first suggested by Poggianti & 
Wu (2000), who also proposed that if the location and thickness of dust depended on the 
age of the embedded stellar populations in a starburst galaxy, then the effects it would have 
on its spectroscopic properties would be quite remarkable. 

Why is it then that only some fraction of dusty starburst galaxies show e(a) spectra? 
Moreover, when observing them in distant clusters, it should not be forgotten that their 
spectrum is an integrated one, comprising a luminosity-weighted 'mean' of the light emitted 
over a large region of the galaxy. Thus while the selective obscuration suggested by 
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Poggianti & Wu (2000) is clearly important, it is the global (100-1000 pc) rather than the 
local (~ 10 pc) processes that need to be investigated in in detail to explain the origin of 
the selective extinction. 

We point out here that a difference in the degree of dust extinction between the central 
(< a few 100 pc) and outer (> Ikpc) regions of a galaxy (~ 1 kpc) could produce the 
above age-dependent extinction effect. To be more specific, if a young stellar population 
is formed in a secondary starburst in the central region where there are greater quantities 
of dust and hence more dust extinction, its contrast in properties to the older and less 
obscured populations in the outer regions of the galaxy will make the extinction appear 
as though it is selective. This idea that the age of the stellar population varies globally 
with radius is not unreasonable and unrealistic, given that secondary star bursts occur 
preferentially in the central parts of galaxies due to the efficient inward transfer of gas 
driven by non-axisymmetric structure (such as stellar bars) and galaxy interaction and 
merging (Barnes & Hernquist 1992; Mihos & Hernquist 1995). Furthermore, numerical 
simulations of the gaseous and stellar distribution in merging disk galaxies (which trigger 
a dusty starburst) have demonstrated that the central young stellar component formed 
by a nuclear starburst, is more heavily obscured by dust than the outer, older component 
initially located in the merger progenitor disks (Bekki, Shioya, & Tanaka 1999). Bekki & 
Shioya (2000) also demonstrated that since younger stellar populations are more centrally 
concentrated and the gas density is higher in the inner region, the mean gas density around 
stars (which is critically important for Ay) depends strongly on stellar age in such a way 
that the density around younger populations is higher than that around older populations. 
Therefore we strongly suggest that galaxies which undergo central dusty starbursts are 
more likely to show selective dust extinction. 

If the suggested radial dependence of stellar age and dust extinction is plausible and 
realistic, are there physical processes operating in distant rich clusters which have an equally 
realistic chance of causing it? Physical processes which are possibilities include: tidal 
galaxy interaction (Icke 1985; Noguchi 1987), tidal effects of cluster cores (Byrd & Valtonen 

1990) , galaxy 'harassment' (Moore et al. 1996), minor merging (Mihos & Hernquist 1994), 
unequal-mass merging (Barnes 1996; Bekki 1998), and major merging (Barnes & Hernquist 

1991) . Given the high probability of high speed galaxy encounters in the rich cluster 
environment, galaxy harassment seems to be the most promising candidate. However it has 
the following three disadvantages in explaining the nature of e(a) galaxies. Firstly, the blue 
cluster galaxies of modest luminosity, i.e. in the luminosity regime where galaxy harassment 
is effective, show predominantly e(b)- rather than e(a)-type spectra (Poggianti et al. 1999). 
Secondly, there is no clear tendency that e(a) galaxies are preferentially located within the 
cores of clusters (Dressier et al. 1999), where galaxy harassment is most likely to trigger 
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a star burst. Thirdly, combined spectroscopic and morphological studies of the galaxies in 
distant clusters (e.g. Couch et al. 1998) have not revealed any significant population of 
starburst objects which has the grossly deformed spiral arm pattern predicted by Moore et 
al. (1996). Of course these observations do not necessarily rule out galaxy harassment as a 
mechanism for forming e(a) galaxies; rather, if galaxy harassment is important, then it will 
produce only low luminosity galaxies located near the cluster core. The second observation 
above would also appear to rule out tidal effects within cluster cores as being important to 
e(a) formation as well. 

A rather high fraction (~ 50 %) of low-2;LlRGs and ULlRGs - considered to be ongoing 
major mergers - are observed to have e(a) spectra (Liu & Kennicutt 1995, Poggianti & 
Wu 2000). The rate of galaxy merging is found to increase very steeply with redshift 
{(X {1 + z) ; van Dokkum et al. 1999) and, furthermore, there is some evidence that this 
trend is seen amongst e(a) galaxies (Dressier et al. 1999; Poggianti et al. 1999). On the 
other hand, the e(a) galaxies found in the distant cluster samples are predominantly disk 
galaxies, which makes it hard to reconcile them with being the products of major mergers, 
which tend to destroy any disk structure. 

We suggest, therefore, that minor and unequal-mass galaxy merging, which is different 
from major merging in that it can leave disk systems intact, is a more likely mechanism 
for reproducing the fundamental properties of e(a) galaxies. Mihos & Hernquist (1994) 
demonstrated that minor merging can excite non-axisymmetric structure in a late-type 
disk, trigger nuclear starbursts, and transform the disk into an early-type system such 
as an Sa or SO galaxy. Furthermore, unequal-mass merging between two disks has been 
demonstrated to create a disk with structure and kinematics strikingly similar to that 
of an SO (Barnes 1996). On that note, Bekki (1998) also found that such merging can 
trigger nuclear starburst, form a very gas-poor disk, and create bulge with characteristics 
remarkably similar to that of typical SOs. Hence as well as contributing to the formation of 
dusty nuclear starbursts and thus to the radial dependence of dust extinction, this more 
minor type of merging may also provide an evolutionary link between late-type disks, c(a) 
galaxies, and early-type disks (SOs) and thus account for the observed decrease in the cluster 
SO fraction with redshift (e.g.. Dressier et al. 1996). An early-type disk galaxy formed 
by a minor/unequal-mass merger, will experience a poststarburst phase (with k-|-a/a-|-k 
spectra), and should contain a bulge remarkable for its young stellar population and thus 
bluer colors (Bekki 1998). Accordingly, a key observational test for assessing the relative 
importance of minor/unequal-mass merging in the formation of e(a) galaxies with selective 
dust extinction, is to investigate the number fraction of disk galaxies with k-|-a/a-|-k spectra 
and prominent blue bulges in distant clusters. 
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4.2. On the origin of less luminous stairburst galeixies with e(b) spectra 

We now focus on the origin of the difference between the e(a) and e(b) spectral classes: 
Why do some starburst galaxies show e(b) spectra while others show e(a) spectra? We have 
demonstrated that the spectral type of a dusty starburst galaxy has a basic dependence 
Avfi in that it will exhibit an e(a) or e(b) spectrum if Avfi is large or small, respectively. 

Poggianti et al. (1999) modeled the time evolution of the global colors and luminosities 
of distant cluster galaxies by using a stellar population synthesis code and found that the 
e(b) galaxies were likely to evolve into sub-luminous systems once the starburst had ceased. 
They also found that the HII region metaUicities of the e(b) galaxies are significantly 
lower than those of galaxies in other spectral classes, and that these low metaUicities were 
consistent with them being low-luminosity galaxies. These observational results suggest 
that one of the physical parameters that can govern the difference in spectral class is 
galactic luminosity. Given our result that the spectral class of a dusty starburst galaxy 
depends on Ayo, this would imply that the luminosity and Ay^Q are somehow related. How 
can this be? One possibility is that for the more luminous galaxies, the thermal feedback 
of supernovae explosions that occur during starbursts is simply not strong enough to blow 
off some fraction of the interstellar gas surrounding the young stars and thus to modify the 
spatial distribution of the gas. The key factor here is the deeper gravitational potential well 
of the galaxy which prevents the gas from escaping the galaxy. As a result of this, the mean 
gaseous density (which is an important factor in determining Avfi) continues to remain 
high, as also its dust extinction. Consequently, the galaxy can show an e(a) spectrum. 

For a less luminous galaxy, most of the dusty interstellar gas can be blown away by 
the thermal feedback effect of the supernovae, as a result of the shallower gravitational 
potential. Consequently, the density of dusty gas around the starburst populations is 
reduced as also the degree of dust extinction. The galaxy thus shows an e(b) spectrum. 
Hence we suggest that owing to the difference in the effectiveness of supernovae feedback 
in galaxies with different luminosities, galactic luminosity can determine the degree of dust 
extinction and optical spectral type. 

Another possibility for the origin of e(b) galaxies is that they are normal star-forming 
systems but with a significantly lower metallicity. Although Poggianti et al. (1999) argued 
that e(b) galaxies must be starburst objects because a spiral-like star formation history was 
not able to produce strong enough [OH] emission, viz. EW([OII])= 40A. However, this was 
based on models with solar metaUicity and indeed Poggianti et al. acknowledged that the 
relation between SFR and [Oil] strength can change with metallicity. In estimating the 
strength of the [Oil] flux for a galaxy with a given star formation rate, the H/3 flux is first 
derived and then used to estimate the [Oil] strength by adopting a certain [0II]/H/3 ratio. 
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Accordingly, if less luminous and less metal-rich galaxies have a larger [0II]/H/3 ratio, then 
they will have larger [Oil] fluxes for a given star formation rate. The large [Oil] fluxes 
observed in e(b) galaxies could therefore be the result of such a metallicity effect rather 
than due to a high star formation rate. Indeed, Kobulnicky et al. (1999) have found that 
the [0II]/H/5 ratio correlates with oxygen abundance (12+log(0/H)) in such a way that 
a galaxy with smaller gaseous metallicity is more likely to show a larger [0II]/H/3 ratio. 
This observational result seems to support the above possibihty that e(b) are just normal 
star-forming galaxies with less metal-rich gas. If this is the case, we do not need to consider 
the reason why there are two different dusty starburst populations. 

4.3. The dusty steirburst nature of the k+a/a+k population? 

Small et al. (1999) found that a significant fraction of the k-|-a population in Abell 851 
{z — 0.42) appeared as sources in deep VLA 1.4 GHz continuum maps of this cluster, and 
suggested that even k-|-a galaxies - classified optically by their lack of [Oil] emission but 
moderately strong RS absorption - can be systems hosting and ongoing dusty starburst. 
Balogh & Morris (2000) also found examples of k-|-a galaxies with Ha emission and, 
accordingly, conchided that some fraction of k-|-a galaxies are likely to be undergoing 
significant star formation. These observations have raised the following two questions: 
(i) Are such k+a galaxies really dusty starburst ones?, and (ii) What determines the spectral 
types of dusty starburst galaxies (or, why are there dusty starburst populations with 
spectral characteristics, viz. e(a) and k+a)? Concerning the first question. Small et al. 
(1999) argued that the radio emission from the k+a galaxies most likely had a starburst 
origin (being the synchrotron radiation emitted from Type II supernovae associated with 
the newly formed population of massive stars) with star formation rates of the order of 
10-30 M0yr~^. Their 'poststar burst' spectral signature at optical wavelengths woTild imply 
that this star formation was completely obscured. Although this interpretation is plausible 
and consistent with other observational results such as color images, submillimeter fiuxes, 
and the [Oil] emission line strengths for these galaxies (Small et al. 1999; Poggianti et al. 
1999), we point out that an alternative interpretation is that there are significant spatial 
variations in star-forming activity across these galaxies, with the central regions in a ' 
poststar burst' phase and the outer regions still in a 'starburst' phase. The galaxy shows 
a k+a spectrum at optical wavelengths because the earlier central star formation (which 
forms the poststarburst population) occurred more intensively than the outer ongoing 
star formation: the integrated spectrum is simply dominated by the light from the central 
poststarburst region. At radio wavelengths, the flux comes exclusively from the outer 
star-forming region. The validity of this interpretation can be assessed observationally. 
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if the detailed two-dimensional spatial distribution of say Ha emission is mapped. 

Furthermore, we can confirm whether Smail et al.'s interpretation of these objects is robust 
by investigating their infrared fiuxes and thereby obtaining more accurate star formation 
rates. Thus it is reasonable to say that future Ho; and infrared observations will resolve this 
first question. 

Concerning the second problem, our present numerical results can provide some clues. 
Our models for starbursts with r = 10^ and 10^ show e(a) and k+a spectra, respectively, 
which implies that the difference in r yields the spectral type difference. The larger value 
of r means that the age-dependence of dust extinction is weaker and thus interstellar dust 
uniformly obscures the starburst populations of different ages. Accordingly, dusty starburst 
galaxies with k+a spectral types are likely to be formed, if younger stellar populations are 
not so preferentially obscured by dust (i.e., the selective dust extinction does not effectively 
work). If the difference in r is really an important factor in explaining the origin of the 
e(a) /k+a dichotomy, the next question is what physical processes drive the r difference. 
Answering this question is obviously beyond the scope of this paper, because r is a free 
parameter of the present one-zone model. However, we suggest that the difference in spatial 
distribution of stars, gas, and dust between dusty starburst galaxies is responsible for the r 
difference. In our future papers, we will investigate how the radial distribution of stars and 
gas affects the spectral types of dusty starburst galaxies. 



4.4. Future work: Spectrodynamical simulations 

Although the present one-zone model of dusty starburst galaxies has clarified some 
essential ingredients of the spectroscopic evolution of distant cluster galaxies, it has the 
following three limitations in analyzing the spectroscopic properties of galaxies. Firstly, 
the degree of dust extinction (represented by Ay), which is one of the most important 
parameters in determining the SEDs of galaxies, is a free parameter in the one-zone 
model. Ay is generally considered to depend on gaseous density, gaseous metallicity, dust 
composition, and the relative distribution of gas and stars (e.g., Mazzei, Xu, & De Zotti 
1992) and so is not calculated in principle by the one-zone model. Secondly, the dependence 
of Ay on the ages of the starburst populations is assumed in an arbitrary way and fixed for 
the one-zone model. The age dependence, which is another important factor in this context, 
depends on the relative spatial distribution of dusty gas and young stars and thus cannot 
be predicted from the one-zone model. Thirdly, the star formation history of a galaxy is 
freely controlled for the one-zone model. The time evolution of the star formation rate is 
suggested to be controlled by dynamical processes such as local gravitational instability 
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in disk galaxies (e.g., Kennicutt 1989, 1998); this is not addressed by the one-zone model 
since it is unable to follow the dynamical evolution of galaxies. In future work, wc suggest 
that the above three properties (e.g., Ay) in the one-zone model should not be freely 
parameterized but predicted from theoretical or numerical studies so that spectroscopic 
properties can be investigated in a more self- consistent manner. 

The most effective way to remove the above three limitations in the classical one-zone 
model of spcctrophotomctric evolution is to derive both the detailed distributions of gas, 
dust, and stars and the time-evolution of the star formation rate, by performing numerical 
simulations of dusty starburst galaxies. Wc have already developed a new code by which we 
can derive photometric evolution of galaxies based on the detailed spatial distribution of gas 
and stars and on the time evolution of star formation rate (Bekki & Shioya 2000). However, 
this code cannot yet enable us to calculate spectroscopic properties of dusty starburst 
galaxies due to the time evolution of gaseous emission not yet being included. Accordingly, 
an important future task is to derive emission line properties in numerical simulations. From 
now on, we refer to these future numerical simulations as 'spectrodynamical simulations', 
by which we can investigate not only spectroscopic properties but also morphological, 
structural, and kinematical ones. These will address the following five questions: (1) Can 
younger starburst populations actually be obscured more heavily by dust than older 
ones, and what physical processes are closely associated with the suggested selective 
dust extinction? (2) What determines the strength of the age-dependence of Ay in dusty 
starburst galaxies and is the strength time-dependent? (3) How do physical processes 
specific to cluster galaxy evolution, such as galaxy harassment (Moore et al. 1994) and 
ram-pressure stripping (Farouki & Shapiro 1984), determine the relative spatial distribution 
of dusty gas and young stars? (4) Are supernovae feedback effects important for explaining 
the spectroscopic properties of less luminous galaxies with e(b) spectra? (5) What is the 
most important dynamical process for the origin of each of the spectral classes observed in 
clusters? 



5. Conclusions 

We have investigated the evolution in spectral properties of dusty starburst galaxies, 
adopting the fundamental assumption that the amount of dust extinction (quantified by 
Ay) depends on the age of the stellar population according to: Ay = Ayx) exp{(T — t)/T} , 
where Ay^ and r are parameters controlling the degree of extinction and T and t represent 
the age of the galaxy and that of its stellar population, respectively. The main conclusions 
from this investigation can be summarized as follows: 
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1. If the young stellar populations (with an age of ~ 10^ yr) are more heavily obscured 
by dust than the older ones (age> 10^ yr) in a dusty starburst galaxy, the integrated 
spectrum of the galaxy is an e(a) type characterized by strong US absorption and 
relatively modest [Oil] emission. The time interval over which the galaxy exhibits 
an e(a) spectrum is fairly long, lasting for ~ 0.7 Gyr of the ~ 1 Gyr duration of the 
starburst. 

2. There arc requisite physical conditions for a dusty starburst galaxy to show an 
e(a) spectrum. For example, the mean Ay^ should be within a certain range 
(~ 1 < Ayfi <~ 10). The age dependence of dust extinction also should be fairly 
strong in a sense that young stars with the ages of ~ 10^ yr are preferentially obscured 
by dust, (i.e., r ~ 10^ yr; more selective). 

3. A dusty starburst galaxy with an c(a) spectrum will evolve into a poststarburst 
galaxy with an a+k (or k+a) spectrum, 0.2 Gyr after the starburst has ended, and 
then into one with a passive k-type spectrum, 1 Gyr after the starburst. This result 
clearly demonstrates an evolutionary link between the following spectral types: e(a), 
a+k, k+a, and k. 

4. The evolutionary path of a dusty starburst galaxy on the EW([0II])-EW(H(5) plane 
depends on the history of star formation within the disk before the starburst, the 
epoch, duration, and strength of the starburst, plus the adopted age-dependence of 
Ay. For example, galaxies of later Hubble type (i.e. with smaller /cdisk) are more likely 
to pass through the bottom-right part of the e(a) region during a starburst. 

5. For a plausible star formation history, a starburst galaxy with the earlier burst epoch 
is more likely to show slightly lower EW([OII]) and slightly higher EW(H5) strengths. 
This result implies that a dusty starburst at higher redshift is likely to be located in 
the bottom-right part of the e(a) region on the EW([0II])-EW(H5) plane. 

6. A dusty starburst galaxy can show an a+k or k+a spectrum even in the starburst 
phase, if the age dependence of Ay is weak (r ~ 10^) for a fixed Ay^. This result 
implies that if young massive stars with different ages (< 10^) are uniformly obscured 
by dust, the galaxy shows an a+k or k+a spectrum. 

7. The duration of a starburst can control the spectral type evolution such that a dusty 
galaxy with a longer starburst time-scale shows a k+a spectrum rather than an e(a) 
spectrum during the starburst phase. This result, combined with the previous point 
implies that the origin of possible starburst k+a galaxies detected by VLA radio 
observations (Small et al. 1999) can result either from the weak age-dependent dust 
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extinction or from the relatively long starburst time scales in these dusty starburst 
galaxies. 

These results clearly demonstrate that the dependence of Ay on ages of young starburst 
populations is a crucially important factor for determining spectral types of galaxies. We 
furthermore stress that the age dependence can be closely associated with the relative 
distribution of dusty gas and young stars for dusty starburst galaxies in distant clusters. 
We will investigate the importance of the relative spatial distribution of gas and stars in 
determining spectroscopic properties of galaxies by using spectro dynamical simulations in 
our future works. 

We are grateful to the referee Michael Dopita for valuable comments, which contribute 
to improve the present paper. Y.S. thanks the Japan Society for Promotion of Science 
(JSPS) Research Fellowships for Young Scientist. K.B. acknowledges the Large Australian 
Research Council (ARC). 
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A. A metallicity dependence of [OII]/H/3 

We used the [0II]/H/9 emission line ratio (adopting a value of 3.01 throughout this 
paper; listed in the table of the PEGASE paper published by Fioc & Rocca-Volmerange 
1997) to calculate the [Oil] emission line strength of a dusty starburst galaxy from the 
derived H/3 line strength in this work. It is important for us to confirm whether or not 
our results, described in the main body of the paper, change if we adopt a [0II]/H/3 ratio 
different to that of the PEGASE model. We therefore demonstrate here that our results do 
not depend strongly on the adopted [0II]/H/3 value, by showing results calculated by the 
photoionization code GLOUDY (Ferland et al. 1998). 

Before presenting the results, we confirm whether our very simple model with the 
GLOUDY code can reproduce reasonably well the [0II]/H/3 values observed for nearby 
galaxies. Figure 11 shows dereddened [0II]/H/9 ratios as a function of oxygen abundance 
(12 + log(0/H) from the observations of Kobulnicky et al. (1999). In this figure, we also 
give our theoretical results based on the GLOUDY code, with an SED of continuous star 
formation of 10 Myr, an electron density of lO^cm"^, and an ionizing parameter U of 10~^, 
for different metallicity models. As is shown in this figure, the GLOUDY-based model can 
successfully reproduce the observational results for galaxies with 12 + log(0/H < 8.5. It 
should be noted here that the [OIIJ/H/3 ratio is less than 3 (corresponding to the value in 
the PEGASE with solar metaUicity) for 8.5 < 12 + logO/H and has a value of ~ 2 around 
solar metallicity (12 + logO/H ~ 8.9). This result suggests that the [Oil] emission line 
strength calculated by PEGASE is appreciably larger (by a factor of ~ 1.5) than that of 
GLOUDY. 

Gonsidering the above difference in the [0II]/H/3 ratio between the two models, we 
investigate the evohition of a dusty starburst galaxy on the EW([0II])-EW(H5) plane for 
Models 1 and 2 by using the GLOUDY code with solar-met allicity (the corresponding 
[OII]/H/? ~ value is 2.08). As is shown in Figure 12, even if we use the GLOUDY code, 
our key results on e(a) formation and spectral evolution do not change (compare Figure 
4 with Figure 12). However, the two models based on the GLOUDY code show lower 
EW([OII]) than the PEGASE-based models during their evolution. This is a result of the 
smaller [0II]/H/3 ratio (2.08), which in turn makes EW([OII]) smaller. As a natural results 
of this, the model with dust shows an e(c) rather than an e(b) spectrum during the actively 
star-forming phase just prior to the starburst. 

For all of the present models, we assumed that the [0II]/H/3 ratio is fixed at a certain 
value for all time steps in a model. This is actually an over-simplification, because the 
gaseous metaUicity (thus [0II]/H/3 ratio) in a galaxy evolves with time as the chemical 
evolution of the galaxy proceeds. Owing to this simplification, the present models show 
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larger values of EW([OII]) in the early evolution phase. In future work, we will consider the 
dependence of the [0II]/H/9 ratio on the gaseous metallicity and thereby construct a more 
realistic and sophisticated model of the spectral evolution of dusty galaxies. 

Lastly, we point out one disadvantage of the present model in deriving the [Oil] line of 
dusty starburst galaxies. Since we used the PEGASE code for the calculation of [0II]/H/3, 
the temperature of the HII region in the present model with the PEGASE code is too low, 
in particular, for high metallicity, owing to the model's not allowing for any dust depletion 
onto grains in the HII regions. As a natural result of this, the [Oil] line strength in the 
present model is severely decreased relative to what it should be. One of the best way to 
avoid this underestimation of [Oil] line strength is to use the [0II]/H/3 derived by Dopita et 
al. (2000) in which dust physics and the depletion of various elements out of the gas phase 
are more properly treated. We have a plan to investigate the [0II]/H/3 of dusty starburst 
galaxies in our future papers by using the model by Dopita et al. (2000). 
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Fig. 1. — The star formation history {upper panel) and the B — V color evolution {lower 
panel) for the fiducial model (Model 1) with the Sb-typc star formation. For comparison, 
the color evolution for the model without dust extinction is also given in the lower panel. 
The starburst begins at Tsb = 7.64 Gyr and stops at Tend = 8.56 Gyr. We assume here that 
the mass of the model galaxy is IO^^Mq (corresponding to the disk mass of the Galaxy). 

Fig. 2. — The spectral energy distribution (SED) with emission hues at T = 1.0 {top), 7.64 
{second from the top), 8.56 {middle), 8.58 {second from bottom), and 13.0Gyr {bottom) for 
the fiducial model. The continuum and the emission lines are represented by the solid and 
dotted lines, respectively. The time T = 7.64, 8.56, and 8.58 represent the beginning of 
the starburst, the end of the starburst, and the epoch when EW([OII]) reaches a minimum, 
respectively. 

Fig. 3. — Upper panel: the evolution of the equivalent width of the [Oil] emission line, 
[EW([OII])], for the fiducial model. Lower panel: The evolution of the equivalent width of 
the absorption line, [EW(H5)]. The solid, dotted, and dashedlines represent the equivalent 
width of the stellar absorption, the equivalent width of emission lines formed in ionized gas, 
and the sum of both lines, respectively. Positive (negative) EW(H5) means that the B.6 line 
is observed in absorption (emission). 

Fig. 4. — The evolution of galaxies on the EW([011])-EW(H5) plane for the fiducial 
model {upper panel) and the model without dust extinction (Model 2). Some points 
(T = 1, 7.64, 8.56, 8.58, 12 Gyr) arc shown as open circles (Model 1) and open triangles 
(Model 2) along the loci of the two models. The classification criteria of Dressier et al. 
(1999) are also superimposed. Both of the models evolve from the e(b) region (T = 1.0 Gyr) 
into the k region {T — 12 Gyr). Note that only the fiducial model can pass through the e(a) 
region during the starburst (7.64 < T < 8.56 Gyr), which implies that the selective dust 
extinction during the starburst is very important for the formation of e(a) galaxies. 

Fig. 5. — The same as Figure 4 but for the models with different Ay^. The results for Model 
3 with Av,o = 1.0, Model 4 with Ay^ = 3.0, Model 1 with Av,o = 5.0, and Model 5 with Ay^ 
— 10.0 are represented by the dotted line with open triangles, short-dashed line with squares, 
solid line with circles, and long-dashed line with pentagons, respectively. The results for the 
five important time steps at T = 1.0, 7.64, 8.56, 8.58, and 13.0 Gyr are indicated by the 
open circles. 

Fig. 6. — As for Figure 5 but for the models with different r. The results for Model 1 with 
T = 1.0 X 10^, Model 6 with r = 3.0 x 10^, and Model 7 with r = 1.0 x 10^ are plotted as a 
solid line with circles, a dotted line with open triangles, and a short-dashed line with squares, 
respectively. 
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Fig. 7. — Top: The time evolution of the star formation rate for models with different 
Tgb- Bottom;. The same as Figure 5 but for the models with different Tsb- The results for 
Model 8 with T^^ = 4.46, Model 9 with Tsb = 6.93, Model 1 with T^b = 7.64, Model 10 
with Tsb = 8.47, and Model 11 with Tgb = 9.45 are represented by a dotted line with open 
triangles, a short-dashed line with squares, a solid line with circles, a long-dashed line with 
pentagons, and a dot-dashed line with crosses, respectively. 

Fig. 8. — As for Figure 7 but for different fc^isk values. The results for Model 12 with 
fcdisk = 0.056 (Sc type). Model 1 with /c^isk = 0.225 (Sb type), and Model 13 with /cdisk = 0.325 
(Sa type), are represented by a dotted line with open triangles, a solid line with circles, and 
a short- dashed line with squares, respectively. 

Fig. 9. — As for Figure 7 but for different ksh values. The results for Model 14 with ksh — 1-0, 
Model 1 with ksb — 0.225, and Model 15 with ksh = 5.0, are represented by a dotted line with 
open triangles, a solid line with circles, and a short-dashed line with squares, respectively. 

Fig. 10. — The same as Figure 7 but for different /final values. The results for Model 16 with 
/final = 0.90, Model 1 with /cgnai = 0.95, and Model 17 with /final = 0.99, are represented by a 
dotted line with open triangles, a solid line with circles, and a short-dashed line with squares, 
respectively. 

Fig. 11. — The observed relation between the [0II]/H/5 ratio and oxygen abundance 
(12 + logO/H) for spiral galaxies {open circles) and metal-poor irregular dwarf galaxies 
(crosses) from Kobulnicky et al. (1999). The solid line represents our numerical results 
based on CLOUDY94 (See the text for more details). 

Fig. 12. — The same as Figure 4 but for the different [0II]/H/5 ratio derived from 
CLOUDY94 (See the text for more details). 
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Table 1. Model parameters 



model 


^disk 


Tsb 


^sb 


/final 


^cnd 


T 


Av,o 






(Gyr) 


(Gyr-^) 




(Gyr) 


(yr) 


(mag) 


1 


0.225 


7.64 


2.25 


0.95 


8.558 


1 X 10*^ 


5 


2 


0.225 


7.64 


2.25 


0.95 


8.558 


1 X 10*^ 





3 


0.225 


7.64 


2.25 


0.95 


8.558 


1 X 10*^ 


1 


4 


0.225 


7.64 


2.25 


0.95 


8.558 


1 X 10^ 


3 


5 


0.225 


7.64 


2.25 


0.95 


8.558 


1 X 10^ 


10 


6 


0.225 


7.64 


2.25 


0.95 


8.558 


3 X 10^ 


5 


7 


0.225 


7.64 


2.25 


0.95 


8.558 


1 X 10^ 


5 


8 


0.225 


4.46 


2.25 


0.95 


5.722 


1 X 10^ 


5 


9 


0.225 


6.93 


2.25 


0.95 


7.920 


1 X 10^ 


5 


10 


0.225 


8.47 


2.25 


0.95 


9.307 


1 X 10*^ 


5 


11 


0.225 


9.45 


2.25 


0.95 


10.193 


1 X 10^ 


5 


12 


0.056 


7.64 


2.25 


0.95 


9.189 


1 X 10^ 


5 


13 


0.325 


7.64 


2.25 


0.95 


8.242 


1 X 10^ 


5 


14 


0.225 


7.64 


1.00 


0.95 


9.870 


1 X 10<^ 


5 


15 


0.225 


7.64 


5.00 


0.95 


8.030 


1 X 10^ 


5 


16 


0.225 


7.64 


2.25 


0.90 


8.148 


1 X 10^ 


5 


17 


0.225 


7.64 


2.25 


0.99 


9.921 


1 X 10*^ 


5 
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EW(H 5 ) (A) 




EW(H 5 ) (A) 




EW(H 5 ) (A) 




EW(H 5 ) (A) 




EW(H 5 ) (A) 




EW(H 5 ) (A) 




EW(H 5 ) (A) 




12 + log (O/H) 




EW(H 5 ) (A) 



